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In this paper, we study a realistic model of quintessential inflation with radiation and matter. 
By the analysis of the dynamical system and numerical work about the evolution of the equation 
of state and cosmic density parameter, we show that this model is a good match for the current 
astronomical observation. The conclusion we obtain is in favor of the model where the modular 
part of complex field plays the role of the inflaton whereas the argument part is the quintessence 
, field. The numerical calculation shows that a heteroclinic orbit (solution of dynamical system) 

interpolates between early-time de Sitter phase (an unstable critical point) and a late-time de Sitter 
' attractor. 
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^ N ■ I. INTRODUCTION 

Cxi ' 

Astronomical observation on the Cosmic Microwave Background(CMB) Anisotropy 1], Supernova Type Ia(SNIa) 
and SLOAN Digital Sky Survey(SDSS) converge on that our Universe is spatially flat, with about seventy 
percent of the total density resulting from dark energy that has an equation of state w < —1/3 and drives the 
accelerating expansion of the Universe which began at a redshift of order one-half. The origin of the dark energy 
remains elusive from the point of view of general relativity and standard particle physics. Several candidates to 
l/" - ) \ represent dark energy have been suggested and confronted with observation: cosmological constant, quintessence with 
a single field £| or with N coupled field Q, phantom field with canonical @ or Born-Infield type Lagrangian Q, 
k-essence @ and generalized Chaplygin gas(GCG) Among these models, the most typical ones are cosmological 
O ■ constant and quintessence which has caught much attention ever since its invention. 

The idea of inflation is legitimately regarded as an great advancement of modern cosmology: it solves the horizon, 
flatness and monopole problem, and it provides a mechanism for the generation of density perturbations needed to seed 
the formation of structures in the universe 01 • In standard inflationary models , the physics lies in the inflation 
potential. The underlying dynamics is simply that of a single scalar field rolling in its potential. This scenario is 
• • ■ generically referred to as chaotic inflation in reference to its choice of initial conditions. This picture is widely favored 
because of its simplicity and has received by far the most attention to date. The properties of inflationary models are 
also tightly constraint by the recent result from the observation. The standard inflationary ACDM model provides 
a good fit to the observed cosmic microwave background (CMB) anisotropy. Peeples and Vilenkin proposed and 
quantitatively analyzed the intriguing idea |12| that a substantial fraction of the present cosmic energy density could 
reside in the vacuum potential energy of the scalar field responsible for inflation (quintessential inflation). After that, 
there were some models to be presented in succession. However, most of these models involve non-renormalizable, 
special potentials and these potentials usually do not have a local minimum, making the conventional reheating process 
unoperative. 

On the other hand, the Peccei-Quinn (PQ) symmetry ^| is the most elegant solution to the stron g CP problem of 
QCD. The global U(1)pq is a spontaneously symmetry breaking (SSB). Weinberg 0] and Wilczek |15| pointed out 
that there is a Nambu-Goldstone boson, "the axion", associated with SSB. Such a global symmetry often arises in 
supersymmetric |16) and superstring-inspired models jl7j . A priori PQ symmetry breaking scale is arbitrary, which 
can be taken a value between 10 2 Gev and 10 19 Gev. In most axion models PQ symmetry breaking happens at a 
temperature T ~ /, when a complex scalar field <&, which carries U{\)pq charge, develops a vacuum expectation 
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value. At temperatures below T <~ / the renormalizable potential for <E> is 



Vx{ri) = -{rf-ff, (1) 

where r\ = ($*<I>)5. However, the argument of ($) is left undecided when the U(1)pq has been broken because of 
Vi(rj) is independent of arg($). The massless £ degree of freedom is the axion: £ = (//A r )arg($), where N is the 
color anomaly of the PQ symmetry ^^|. The axion is massless at T ^ QCD scale, but the axion develops a mass due 
to instanton effects at low temperatures. In the N = 1 case, the axion mass is easy to concretize by using a potential 



U 2 (£) = Af 4 [l-cos(£//)], (2) 

where M 4 = ?n|/ 2 and the axion develops a mass mt. The conception that the modular part of a complex field with 
PQ symmetry may drive inflation is also not new jlflj . Another version of a complex field was proposed as a realistic 
candidate for inflation l20j and, separately, for quintessence |2ll in a natural fashion. The two-axion model has also 
been suggested for quintessence [23] or inflation 23] . The axion model can also be used for quintessence |24} and, 
separately, for inflation |2"j| from extra dimensions. 

Recently, a fascinating model for quintessential inflation has been proposed by Rosenfeld and Frieman (RF) j2(J 
with aid of the axion theory. In their model, a renormalizable complex scalar field $ described by a lagrangian with a 
U(l) symmetry spontaneously broken at a high energy scale / ~ 10 19 Gev and explicitly broken by instanton effects 
at a much lower energy can account for both the early inflationary phase and the recent accelerated expansion of the 
Universe. The modular and argument parts of the field <& were identified with the inflaton and the quintessence fields 
respectively. 

In this paper, we put the emphasis on the study of the dynamics of RF model with radiation and matter composition. 
According to the phase space analysis and numerical calculation, we show that how the modular part of complex 
field plays the role of the inflaton whereas the argument part is the quintessence field. Phase space methods are 
particularly useful when the equations of motion are hard to solve analytically for the presence of radiation and dust 
matter. In fact, the numerical solutions with random initial conditions are not a satisfying alternative because of these 
may not reveal all the important properties. Therefore, combining the information from the critical points analysis 
with numerical solutions, one is able to give the complete classification of solutions according to their early-time 
and late-time behavior. There are stable and unstable critical points for RF model with radiation and dust matter 
composition. Critical points are always exact constant solution in the context of autonomous dynamical systems. 
These points are often the extreme points of the orbits and therefore describe the asymptotic behavior. In our case, 
the unstable critical points are corresponding to a de Sitter behavior at very first moments of the Universe. The 
stable critical points are corresponding to the second stage of accelerated expansion began at a redshift z ~ 1.5 and is 
still operative, which alleviates the fine tuning problem. The equation of state w varies with the cosmic evolution and 
approaches towards -1 asymptotically showing the existence of a cosmological constant at late times. If the solutions 
interpolate between critical points they can be divided into heteroclinic orbit and homoclinic orbit. In our case, the 
heteroclinic orbit connects unstable and stable critical points. The numerical calculation shows that brought the 
Universe back to the usual Friedman-Robertson- Walker expansion, then the second stage of accelerated expansion 
began at z ~ 1.5. 



II. PHASE SPACE OF U(l) QUINTESSENTIAL INFLATION 



For convenience, we investigate firstly the global structure of the dynamical system without radiation and matter. 
The equations of motion for RF model [2(j are 

rj + 3Hf l -^r 1 + V{(T 1 )=0, 

' . 2 (3) 



and the Einstein equation is 



H 2 = y(p„ + Pc), 



(4) 
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where H equals a/ a is Hubble parameter and n 2 — 8irG. Here we define the energy densities: 

To gain more insights into the dynamical system, we introduce the new variables 

u = 77, x — £, 
v = ri, y = i, 

then the equations of motion could be reduced to 

du 

dv 
~d~b 
dx 
~dl 
dy 
dt 



(5) 



(6) 



-3Hv+^--V{(u), 



(7) 



y- 



2v f 2 
-(3H+-)y-^V 2 \x) 



And H could be rewritten as 



H 



K 1 



1 u 



., yi v 2 + V 1 (u) + --y 2 + V 2 (x) 



1/2 



(8) 



Here we have V x (u) = ±(u 2 - f 2 ) 2 and V 2 (x) = M 4 [l - cos(x//)]. 

Before we carry out the numerical study, we would like to analyze the system of equations qualitatively. The 

critical point of the autonomous system Eqs.lQ) is (u c , 0, x c , 0), where defined by V{(u c ) = 0,V^(x c ) = 
respectively. So one can obtain the critical points: 



(u c , v c , x c , y c ) = (±f, 0, kfir, 0), 



(9) 



where k is any integer. 

In order to investigate the stability of these critical points, we can linearize the Eqs.Q about the critical points 

y = Ay, (10) 
df 

A = -w=(u c , v e ,x c ,y c ). (11) 
dy 

Here y = (it, v, x, y) T and / is short for the right hand sides of Eqs.Q. Then we have the stability matrix A 



/ _J 

-V{'(u c ) -y/3K 2 V 2 (x c ) 


\ 

The eigenvalues of the stability matrix are 



r 1 

2« 






° \ 











1 


-V>'{x c ) 


-^3k 2 V 2 {x c )J 


3k 2 V 2 (x c ) 


-4^"K), 


3k 2 V 2 (x c ) 


-4^"K), 


3k 2 V 2 (x c ) 


-4F 2 "(x c ), 



(12) 



-(-^3k 2 V 2 (x c ) + yj3K 2 V 2 (x c ) - W 2 "{x c ) 
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It is clear that the critical points are always stable when k is even because the real part of all the eigenvalues is 
negative. So these critical points correspond to late-time attractor solutions. When k is odd, the fourth eigenvalue 
in Ea. ()12jl is always larger than zero and the corresponding critical points are unstable, which are exact constant 
solution with de Sitter behavior at early-time of the Universe. This is to say that the dynamical system has stable 
critical points at the minimum of the potential Vx(rj) and the minimum of the potential V2(£). Next, let's read out 
the physical implications when the system is at the critical point regime. The cosmic parameters are 



0, = ^, n £ = £ (is) 

Pc Pc 



and the equation of state of the complex field is 



v 



2 + ^y 2 -2V 1 (u)-2V 2 {x) 



w = 7T, ; ■ (14) 

Clearly, from Eas. (|13(l and flj 1 4|> . one can find that w = — 1, Q, n = and f2j? = 1 at the late time attractor. 
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FIG. 1: Evolution of the equation of state at the presence of radiation and matter. The curves correspond to the different initial 
conditions. 




FIG. 2: Evolution of cosmic density parameter of quintessence energy density Q^, inflaton energy density fl v , radiation energy 
density f2 r and matter energy density Q m .The curves correspond to the different initial conditions 
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III. EVOLUTION AT THE PRESENCE OF RADIATION AND MATTER 



If wc arc to search for a realistic model, we must know exactly how the complex scalar field <& couples to ordinary 
matter. For convenience, we assume that coupling is very weak following traditional ideas. In fact, the most significant 
feature of all the axion couplings is that they are proportional to N / f: the larger the PQ SSB scale, the weaker the 
axion couples. Furthermore, we have taken the / ~ 10 19 Gev so that couples can be neglected. Certainly, the ordinary 
matter still affect the evolution of <E> via their contribution to the general expansion of the Universe. 

Now we solve the equations of motion via a numerical calculation at the presence of radiation and matter and 
obtain the results that will confirm our qualitative analysis. To do this, we introduce the new variable 



Therefore, the Eqs.JJJ) are reduced to 



N = \na. (15) 



du v 

In ~h' 

dv_ ^ uy*_ _ V{(u) 

dN + pH H ' 

(16) 

dx y 
dN = H' 



dN~ {6+ uH )V ■ 



and 



H 2 = H?E 2 (N), (17) 
where Hi denote the Hubble parameter at an initial time. E(N) is defined [23] as 



£(V,= .M* {^ v2 + Vl{u) + \$ y2 + V2{x) 



1/2 



(18) 



Q rt i and VL m ^ are the cosmic density parameters for radiation and dust matter at the initial time. Note that when N 
goes to be very large, or at late-time, the contribution to H from matter and radiation will become negligible. 

Solving these equations will give us some insights into the evolution of the fields and the quantities of interest. The 
corresponding conclusions of numerical calculation are shown in Fig^and FigEl This numerical solution describes a 
heteroclinic orbit, which interpolate between an unstable de Sitter critical point and a late-time de Sitter attractor. 
Therefore, this model is more natural to explain the two stages of acceleration. A point worth emphasizing is that 
the ordinary matter (radiation and dust) affect the evolution of the scalar field via their contribution to the general 
expansion of the Universe because of the couples can be neglected in the model. In Fig^ the behavior of the equation 
of state parameter is shown. At initial time, w = — 1 corresponds the de Sitter expansion (inflationary phase), then 
it increases and becomes positive. After arriving at the value 1/3, the Universe comes to the radiation dominated 
epoch and w stays on a broad platform. Next, w drops to zero and stays on a narrow platform. Finally, w drops 
below zero and approaches to —1, which corresponds second stage of accelerated expansion. The evolution of cosmic 
density parameters are shown in FigEI The argument part contribution stays at f2 = at the very first moments 
of the Universe and becomes 1 in late-time, which plays the role of quintessence field. On the other hand, the modular 
part contribution Q v plays the role of inflaton. During the whole evolution, the radiation energy density and the dust 
matter energy density become dominate respectively. Therefore, we see that the constraints arising from cosmological 
nucleosynthesis and structure formation are satisfied. 



IV. CONCLUSION AND DISCUSSION 



In this paper, we have discussed the cosmological implication of a complex scalar model of quintessence inflation 
with radiation and matter. Using the numerical calculation for this model, we show that the argument contribution 
is negligible at the early epoch of the universe while it becomes dominate with the time evolving and the evolution 
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of the modular part is the opposite, which is a viable way to unify the description of the inflationary stage and the 
current accelerated expansion. Analysis to the dynamical evolution of the complex scalar model indicates that it 
admits a late-time attractor solution, at which the field behaves as a cosmological constant. Obviously, attractor and 
heteroclinic orbit are both alleviate the fine tuning problem. We can also extend our discussion to the lagrangian 
with non-abelian symmetry, which is an interesting extended work and we will discuss in another preparing work. 
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